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1 Introduction 

We consider the following first order systems of mathematical physics. 

1. The Dirac equation with scalar potential. 

2. The Dirac equation with electric potential. 

3. The Dirac equation with pseudoscalar potential. 

4. The system describing non-linear force free magnetic fields or Beltrami 
fields with nonconstant proportionality factor. 

5. The Maxwell equations for slowly changing media. 

6. The static Maxwell system. 
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We show that all this variety of first order systems reduces to the equa- 
tion 

Df + f-S = 0, (1) 

where D is the Moisil-Theodoresco operator (introduced by Hamilton) act- 
ing on biquaternion valued functions / according to the rule Df = J2k=i e kdkf, 
3^ = e k are standard quaternionic imaginary units; the function / of 

real variables x\, x 2 , x 3 has the form / = Ylk=o fk e k, fk S C, k = 0, 1, 2, 3 
and a is a purely vectorial biquaternion valued function. We reduce the 
solution of equation Q to the solution of a Schrodinger equation with bi- 
quaternionic potential. In some important situations the biquaternionic 
potential can be diagonalized and converted into scalar potentials. 

2 Notations 

We will consider the algebra H(C) of complex quaternions or biquaternions 
which have the form q = J2k=o l^k where {q k } C C, eo is the unit and 
{e^l k = 1, 2, 3} are the quaternionic imaginary units, that is the standard 
basis elements possessing the following properties: 

e o = e o = -ejfcj e Q e k = e k e = e k , k = 1, 2, 3; 

ei e 2 = -e 2 e 1 = e 3 ; e 2 e 3 = -e 3 e 2 = e.\\ e 3 e 1 = -eie 3 = e 2 . 

We denote the imaginary unit in C by i as usual. By definition i commutes 
with e k , k = 0, 1, 2, 3. 

The vectorial representation of a complex quaternion will be used. Namely, 
each complex quaternion q is a sum of a scalar q and of a vector q : 

q = Sc(g) + Vec(g) = <? + ~<f , 

where "g* = Ylk=l 1k e k- The purely vectorial complex quaternions (Sc(g) = 
0) are identified with vectors from C 3 . The following conjugation operations 
will be needed. The quaternionic conjugation is defined as follows q = qo—~q, 
the complex conjugation: q* = Reg — ilmq and the following involutive 
operation q^ = e k qek which changes signs of two components, for example, 
<7 (1) = -eiqei = q e + qiei - g 2 e 2 - q 3 e 3 . 

By M p we denote the operator of multiplication by a complex quaternion 
p from the right-hand side: M v q = q ■ p. 

We will intensively use the fact that the algebra of complex quaternions 
contains a subset of zero divisors & which are characterized by the equality 
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go = ~Q 2 , where ~~q 2 = — ( q , q ), or equivalently g 2 = 2gog. Hence if g 6 6 
and go = 1/2 then g is an idempotent. Using this fact we introduce the 
following multiplication operators P^ = ^M( 1=tiefe ), k = 1, 2, 3. It is easy to 
see that for each k the operators and P k ~ represent a pair of mutually 
complementary, orthogonal projection operators on the set of EI(C)-valued 
functions. More information on the structure of the algebra of complex 
quaternions can be found for example in ^Jj or 

Let / be a complex quaternion valued differ entiable function of x = 
(xi, X2, X3). Denote 

k=i K 

This expression can be rewritten in vector form as follows 

Df = - div / + grad / + rot / . 

That is, Sc(Df) = — div / and Vec(Df) = grad/o + rot / . Let us notice 
that D 2 = — A. The operator D + M a we will denote also by D-$. 

Let us introduce an auxiliary notation / := f(x\,X2, —xz). The domain 
Q, is assumed to be obtained from the domain 1] C I 3 by the reflection 

3 First order systems reducing to equation (0Q) 

3.1 The Dirac equation with scalar potential 

Let V denote the classic Dirac operator for a free particle with a specified 
energy w£l 

3 

V = iu^o + ^2 j k d k + im. 
k=l 

Here 7^, j = 0, 1, 2, 3 are usual 7-matrices (see, e.g., 0) and m G R. 
The Dirac operator with scalar potential has the form (see, e.g., |23| ) 

V sc = V + i(p sc I, 

where <p sc is a scalar real-valued function of x and / denotes the identity 
operator. 

In |13j (see also |17j . |19j ^ the following transformation was introduced. 
A function $ : ft C M 3 — > C 4 is transformed into a function F : £1 C M 3 — > 
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H(C) by the rule 

F = A[<$>] := - (-($1 - 5> 2 )e + i(<S>o - <& 3 )ei - ($ + $3)e 2 + + $2)e 3 
The inverse transformation A^ 1 is defined as follows 

$ = A-\F] = {-iFx - F 2 , -F - iF 3 , F - iF^m - F 2 ). 

Let us present the introduced transformations in a more explicit matrix form 
which relates the components of a C 4 -valued function <3? with the components 
of an H(C)-valued function F: 



F = A[&] = 

and 

<f> = A- 1 [F] -- 
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The following equality is valid 

D^ sc = ~A 7m V sc A-\ (2) 

where D-^ sc = D + M~^ sc and 7x sc = —(iuei + (m + (p sc )e 2 ). Thus the 
function is a solution of the Dirac equation with scalar potential 

V SC <S> = in n 

if and only if the function F = A& is a solution of the equation 

D-g F = in fi. 

3.2 The Dirac equation with electric potential 

The Dirac equation with electric potential has the form 

V d = V + 1^,70, 
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where cp e i is a real- valued function. 
We have an equality similar to ©: 

D^ el = -^ 7 i7273^ e ^- 1 , (3) 

where -D-^ eI = D + M° el and ~a e \ = — (i(ui + </?ez) e i + me 2)- Thus the 
equation V el Q = is equivalent to the equation D-g F = 0. 

3.3 The Dirac equation with pseudoscalar potential 

The Dirac equation with pseudoscalar potential has the form (see, e.g., [23] ) 

VP S = V + <^ s7 o75, 

where (p ps is a real-valued function. We have the following equality (see 
[T7|). similar to © and ©: 

D + vl + = -^7i727 3 2> ps >l _1 , 

where v = —i(p ps and = —{iuje\ + me?). 

Suppose that /? ^ 6, that is m? ^ w 2 and denote S ± = -^M^ x±l3 \ 
where the complex number A is chosen such that A 2 = f3 2 . We have |14j 

D + vl + = S + (D + + A)/) +S~(D + (y- A)/). (4) 

The operators of multiplication S + and S~ are mutually complementary 
projection operators on the set of H(C)-valued functions and commute with 
the operators in parentheses in (|4"|). Thus / satisfies the equation 

(D + J/I + M~&)f = (5) 

if and only if the functions / + = S + f and /~ = S~f are solutions of the 
equations 

{D + u + X)f+ = (6) 

and 

(D + v — X)f~ = (7) 

respectively. In other words, given /+ and /~ solutions of © and (J7J), the 
function / = S + f + + S~ f~ will be a solution of ©. 
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Equations (jSJ) and (JJJ) have a quite convenient form for studying, nev- 
ertheless in order to reduce them to the universal form Q we make one 
additional step. We have 

D + (y + A)/ = P+ (D + M {u+X ^ ) + Pf (D - M iu+ ^ iei ) (8) 

and 

D + (y - A)/ = P+ (D + Af(" _ A)tei ) + Pf (-D - M^'^ 1 ) . (9) 

Thus / + and / _ are solutions of © and (J7J) respectively if and only if the 
functions f ++ = P^~ f + and / h = Pf f + are solutions of the equations 

(D + M ( - U+X ^ iei )f ++ = 0, (10) 



(D-M^ u+X ^ ei )r + = 0, (11) 

and the functions Z" 1 = P^f~ and = Pf /~ are solutions of the 

equations 

(D + M^^ 61 )/^ =0, (12) 
(D - M^-V' 161 )/— =0. 
The obtained result we resume in the following statement. 

Proposition 1 For (3 ^ & and A 2 = f3 2 the following equality is valid 

D + vI + M? = P+S + {D + M^ +A ) iei ) + Pf S + (D - M^ +A ) iei ) 
+ P+S~(D + M^ A ) iei ) + Pf S~(P> - M^- A ) iei ) 

which implies 

ker(D + vl + M~@) = P+5+ ker(D + M^ +A ^ ei ) Pf S+ ker(D - M^ +A ) iei 
© ker(L> + M^~^ iei ) © Pf5~ ker(D - M^~^ iei 

where ker means the set of null- solutions in a domain of interest Hcl 3 . 

In this way equation (JSJ and hence the Dirac equation with pseudoscalar 
potential reduce to four equations of the form (jT|). 
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3.4 Force-free magnetic fields 

Force- free magnetic fields appear as an important class of special solutions of 
nonlinear equations of magnetohydro dynamics and are intensively studied 
in different branches of modern physics (see, e.g., PQ, [H], [S]> E3) [TTj . 
US, HQ, HH, QUI > US])- They are characterized by the following pair of 
equations 

divB = (13) 

and 

rotB + ^B = 0, (14) 
where v is a scalar function. This system can be rewritten in the form 

(L> + z;)B = 0. (15) 

It will be more convenient for us to extend the class of solutions of Q15j) and 
to consider not only its purely vectorial solutions but more general complete 
biquaternionic functions. Thus we consider the equation 

(D + u)f = 0, (16) 

where / is an HI(C)-valued function. Solutions of (|13|) . (|14|) represent a 
subset of solutions of ()16|) which fulfil the additional requirement Sc/ = 0. 

Now by analogy with the preceding subsection we obtain that / is a 
solution of (|16|) if and only if the functions / + = f and /~ = P{~ f are 
solutions of the equations (D + M iuei )f + = and [D - M iuei )f~ = 
respectively. 

Thus the system ((T3*|) . (fTl|) reduces to a pair of equations of the form 

©• 

3.5 Maxwell's equations for slowly changing media 

Consider first the general Maxwell system 

rotH = ed t E + j, (17) 



rotE = -^H, (18) 
div(eE) = p, (19) 
div (pH) = 0, (20) 
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where e and \i are assumed to be functions of spatial coordinates only. It 
can be rewritten in the following form 

(D + M t )E = --d t H--£ = , (21) 
c ye 

and 

(D + M^)H = - c d t E + v77j, (22) 

where £ = y^E, H = yTIH, c = 1/y^/Z, = and 7? = ^^E. 

In a sourceless time-harmonic situation we obtain the equations 

D-gE = ivH and D-^H = -ivE . (23) 

Here v = u/c. 

The medium is said to be slowly changing when its properties change 
appreciably over distances much greater than the wavelength [21 El]- Usually 
this is associated with the possibility of reducing the Maxwell equations (|2*3*)l 
to the Helmholtz equations 

(A + v 2 )E = and (A + u 2 )H = 0. 

It is easy to check that such a reduction is possible if only \~e\ and | \x | are 
considered as relatively very small and the terms containing the vectors e 
and fx are supposed to be negligible. Then (|23|) take the form 

DE = ivH and DH = -ivE 

and can be diagonalized. For the functions Tp = E + iH and ip = E — iH 
we obtain the equations 

(D - v)~$ = and (D + v)~^} = 0. 

Now, by analogy with equation (|T5|) each of these two equations can be 
rewritten in the form [|T|). 

3.6 The static Maxwell system 

When the vectors of the electromagnetic field do not depend on time, from 
1)21(1 and (|2*2*)l we obtain two independent equations 

(D + M^)E =— (24) 
v e 

and 

{D + M^)H = ^Jiy (25) 
In a sourceless situation both equations reduce to 
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3.7 Two important cases 

Let us resume the results presented in this section. We see that the consid- 
ered here six physical models reduce to equation (^Q). In the case of the static 
Maxwell system a is a gradient of some scalar function. In all other cases 
in general a must not be necessarily a gradient, though such a possibility is 
not excluded. In the first five models only one of the components of a is a 
function while the other two are constants. Thus we are interested basically 
in the following two situations. 

1. ~a is a gradient of some scalar function ip: ~~a = X7ip. 

2. ~a has the form ~u = a.\{x\, X2, xs)ei + 02^2 + a^e^, where ot\ is a 
complex valued function and 02, 03 are complex constants (of course, 
when qi = a\{x\) we have the first case again). 

4 Factorization of the Schrodinger operator 
4.1 Scalar Schrodinger equations 

Consider the Schrodinger operator — A + v applied to a scalar function <j). 
Here v is some complex valued function. Let ~a be a purely vectorial bi- 
quaternion valued function such that 

D~a + (a f = -v. (26) 

Then as was observed in jS] and 0], the following equality is valid 

(-A + v)<(> = (D + M^)(D - M^)^. (27) 

This equality gives a certain relation between solutions of (0) an d null- 
solutions of the Schrodinger operator. Moreover, a fundamental solution of 
the operator D-g can be obtained if a fundamental solution of the Schrodinger 
operator is given. In some cases, for instance when ~a = ai(xi)ei, this leads 
to construction of integral representations for solutions of (^Q) (see [T5]). 

Equation (|26|) represents a generalization of the famous Riccati differ- 
ential equation. In ^] (see also ^7]) it was shown that it is not merely a 
formal resemblance. On the contrary, the well known Euler theorems, the 
Weyr theorem as well as the Picard theorem descibing the unique properties 
of the Riccati equation were generalized for equation (|26|) . 

Equation ()26|) necessarily implies that of is a gradient, because from the 
vector part of ()26|) we have that rot ~ct = 0. Thus it can be useful only for 
the first situation from Subsection 3.7. 
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Let us consider the product of the operators D-g and D_-$ in appli- 
cation to a twice differentiable biquaternion valued function u and for any 
differentiable purely vectorial biquaternion a . We have 

3 

D^D_^u={D + M^){D-M^)J2^kek 

k=0 

3 

= J2( Mek ( D + M« W )(D - M" (fe) )ufe)- 
fc=o 

For each component Uk we use (|27j) and note that (~a = ~a 2 . We obtain 

3 

D^D^u = ^(-Au fc - ^ 2 u k - {D^)u k )e k . (28) 

fc=0 

From (|2"8|) we see that the equation D-gD_-gu = is equivalent to four 
scalar Schrodinger equations if only D~a^ e ' is scalar for any k = 0, 1, 2, 3. It 
is easy to verify that it is possible if only a has the following form 

~a = ai(xi)ei + a 2 (x 2 )e 2 + a 3 (x 3 )e3. (29) 

Let us consider this case in detail. First of all we notice that for the Dirac 
operator with scalar, electric or pseudoscalar potential this restriction on a 
implies that the potential is an arbitrary function of one spatial coordinate. 
In the case of force-free magnetic fields (Subsection 3.4) the proportionality 
factor v is a function of one variable. In the case of the electromagnetic 
field in a slowly changing medium (Subsection 3.5) the wave number v is a 
function of one variable. Such media are known as stratified media. 

For the static Maxwell system div(eE) = 0, and rot E = we obtain that 
it is equivalent to Q with a having the form (|29j) iff the permittivity e has 
the following form e = £i(xi)e2(^2)e3(^3)) where e±, e 2 and £3 are arbitrary 
functions, otk are related to £k in the following way oik = (dkSk)/(2ek), 
k = 1,2,3. 

Thus for all considered physical models ~a in the form (|29[) corresponds 
to quite interesting, nontrivial situations. Let us study in detail the structure 
of solutions of equation Q when a has the form (|29j) . 

Denote 

v k = -DTtW - ~a 2 , (30) 

and 

w k = D~a (fc) - "a 2 , fc = 0,l,2,3. (31) 
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Proposition 2 Let f be a solution of £7J) with a having the form \2y\) . 
Then the components f k are solutions of the Schrddinger equations 

(-A + w k )f k = 0, k = 0,1,2,3. (32) 

Proof. Assume that / is a solution of (pQ) . Then considering the equation 
D_-^D-gf = we arrive at the equations (|3*2*)) for the components fk- ■ 

The following statement gives us a method for constructing exact so- 
lutions of (^Q) having obtained solutions of the corresponding Schrddinger 
equations. 

Proposition 3 Let a be of the form and four scalar functions gk, 
k = 0, 1, 2, 3 satisfy the following equations 

(-A + v k )g k = 0. (33) 

Then the function 

f = (D-M^)g (34) 
is a solution of Q), where g = Ylk=o9k e k- 

Proof. This is an immediate consequence of (|28|). ■ 

Example 4 Let ~a = Yl\=i Xk -b k &k > w ^ iere °k ar & arbitrary complex con- 
stants. We have 

Vo = -D~a - ~a 2 = 0. (35) 

This means that taking any scalar harmonic function go we are able to con- 
struct a solution of 0) in the form f = {D — M a )g$. Calculating v\ , V2 
and V3 according to \3U\) we obtain 

(36) 
(37) 

v*= 2 ( tz 1 -TY2 + 7Z (38) 



(xi - 61) 2 (x 2 - b 
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We will not try to find general solutions of the Schrodinger equations 
with these potentials. We show instead how one can always obtain a class of 
exact solutions of (J2) and of the corresponding Schrodinger equations P3*|) 
when a has the form Q29JI . 

Let us look for a one-component solution of ©: / = f k e k . We have 
Df k e k + f k aWe k = 0. That is 

(39) 

Jk 

It is easy to see that for a of the form (|29j) the functions f k are 

j _ e -(Ai+A2+A 3 ) t _ e -Ai+A 2 +A 3 j 2 _ e Ai-A 2 +A 3 t _ e Ai+A 2 -A 3 

where A k = A k (x k ), k = 1, 2, 3 is an antiderivative of a k . Thus the function 
Ylt=o c kfk&k is a solution of Q, where c k are arbitrary complex constants. 

Concerning the Schrodinger equations (|33|) it is useful to consider the 
quaternionic Riccati equation (|26|) . Any solution of it can be represented in 
the form ~a = — ^, where tp is a solution of the equation — A</3 + u<£ = (see 

H2|, H3). For "a from (j^ we obtain that ~a = where (^ fc = l/f k 
and consequently the functions tp k are solutions of (|33^ for a corresponding 
fc. 

Let us see what are the functions f k and (p k in Example |IJ 
Example 5 For a from the preceding example we have 

^ = 1 . = (x 2 - b 2 )(x 3 - b 3 ) 

{x\ - h)(x 2 - b 2 )(x 3 - 6 3 )' xi-bi 

(xt - bx)(x 3 - 63) , _ {x\ - bi){x 2 - b 2 ) 

J 2 — 7 ) J 3 



x 2 -b 2 x 3 - b 3 

^3 

x\ - bi 



Hence the function f = Y%=n c kfk e k is a solution of 0) and the functions 



ip = (xx - bi){x 2 - b 2 )(x 3 - 63), ipi 



(x 2 - b 2 )(x 3 - 63) ' 



x 2 -b 2 x 3 - b 3 

^2 = 7Z — rT7z — rr> w 



(xi - h)(x 3 - 63) (xi - h)(x 2 -b 2 ) 

are solutions of h'3'J\) with potentials h'35\) - H33\) respectively. 

Proposition 6 Let ~a be of the form \29\) and f = (D — M a )g a solution 
of 0), where g = Yl\=Q9k^k- Then necessarily g k are solutions of ]3ty 
respectively. 
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Proof. Consider 

3 

(D + M% = (D + M«)(D - M°)^2g k e k 

k=0 

3 

= ^2(M ek {D + M* W )(D - M" (fc) ) 5fc ) 

k=Q 
3 

= J2(M e k(-A + v k )g k ). 

k=0 

From assumption of the proposition we obtain that (— A + v k )g k , k = 
0,1,2,3. ■ 

Let Q, be a domain in R 3 which in particular may coincide with the whole 
space. Let F(£l) and G(£l) be some functional spaces. The lH(C)-valued 
function / is said to belong to a functional space if each of its components 
f k belongs to it. 

Proposition 7 Let~a be of the form \2y\) . Assume that the equation (— A+ 
w k (x))u(x) = ^(x), xe!l, k = 0, 1,2, 3 is solvable for any /j, G F(£l) and 
the solution u belongs to G(Q). Then the equation 

(D - M 1} )g = f (40) 

is solvable for any f G F(Q) and the solution g belongs to \mD-^{G{Vt)) . 

Proof. Let u k , k = 0, 1, 2, 3 be solutions of the equations (— A+w k )u k = 
f k . Then g = (D + M") Yll=o u k e k is a solution of (|UJ). ■ 

Remark 8 We do not specify the functional spaces here because the results 
on solvability of the inhomogeneous Schrddinger equation are numerous and 
correspond to very different situations. Let us give an example. 

Example 9 Let Q, = M? and w has the form 

w = w — m 2 , (41) 

where w G C^°(R 3 ), m > 0. The equation (—A + w)u = (j, is uniquely solv- 
able 0/ for any /u G L 2 , a and u G #f oc (R 3 ), u = O(^), ^ - imu = o{ ± ) 
for |x| — > oo. Here L 2ya denotes the space of square integrable functions with 
support in a ball of radius a. 

Thus from Proposition ^ we have the solvability of \4-Uj ) for any ~a 
such that w k defined by \31\) have the form \41\j , for example, for ~a = 
ELi Pk(xk) + im k , where f3 k G C^°(IR 3 ) and m k > 0. 
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Proposition 10 Under the conditions of Proposition^ any solution of 0) 
from F(Q) has the form 

f = (D- M")g, (42) 
where g = Ylk=o9k e k- an d 9k satisfy the equations in Q. 

Proof. From Proposition [7| it follows that / can be represented in the 
form ([42 j). where g G im D-g(G($l)). From Proposition El we obtain that 
are solutions of PH)), ■ 

For the Schrodinger operator there are developed different approaches for 
obtaining asymptotic fundamental solutions under some additional assump- 
tions. A fundamental solution can be used for construction of a right-inverse 
operator for the Schrodinger operator which gives a possibility to solve in- 
homogeneous equations. In the following proposition we show that having 
constructed right-inverse operators for the Schrodinger operators —A + v k , 
k = 0, 1, 2, 3, one can construct a right-inverse operator for D-g. 

Proposition 11 Let a be of the form \29^ and T k such operators that 
for any ip £ F(£l): (—A + v k )T k ip = <p in £1, k = 0,1,2,3. Then for 
any f = J2k=o fk e k G F(tt) we have D^T-gf = f in O, where T-gf = 

(D-M«) (zl=o(Tkfk)e k ). 
Proof. Consider 

3 

(D + M")T-/ = ^(M efc (£> + M 7l(k) )(D 

k=0 

3 

= ^2(M^(-A + v k )T k f k ) 

k=0 



4.2 Schrodinger equations with quaternionic potentials 

Let us assume that a has the form 

"a = ai(xi,X2,x 3 )ei + a 2 e 2 + oj 3 e3 (43) 

where a\ is an arbitrary complex valued differentiable function and ot2, «3 
are complex constants. In Section 3 we saw that the first five physical models 
considered here correspond to this situation. 



- M^ k) ){T k f k )) 

3 

= fk e k = /• 
k=0 
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Note that in this case D~~a = D~a^ = —Wa^ = —D7y^ 3 \ Hence from 
(|28[) we obtain 

D-^D_-gu = -Au - ~a 2 u - (D~a )(ii e + - u 2 e 2 - u 3 e 3 ), 

where u = Efc=o u fc e fc- 

Denote Cu = —e\ue\ = u^e® + u\e\ — u 2 e 2 — u^e^, Au = —Au — ~a 2 u 
and Bu = — {D~a)u. That is 

DuD-S = A + BC (44) 

and 

D_^D^ =A-BC. 

Denote Q ± = \{I±ie 1 C). 

Proposition 12 Solutions of equation 0) (with ~a of the form \4^ ) have 
the form f = D_-^u, where u = Q + v + Q~w and v, w are solutions of the 
following Schrddinger equations with quaternionic potentials 

(A + Biex)v = (45) 

and 

(A - Bie^w = (46) 

respectively. 

Proof. First, we notice that Q ± B = BQ ± and Q ± C = \[C ± ie x ) = 
\{Cie\ ± I)ie\ = iQ^iei. 

Applying Q + and Q~ to the equation 

(A + BC)u = (47) 

we see that it is equivalent to the pair of equations (A + Bie±)Q + u = 0, 
{A— Bie\)Q~ u = 0, and Q + , Q~ commute with the operators in parentheses 
and represent a pair of mutually complementary projection operators on the 
space of EI(C)-valued functions. Thus we have that u is a solution of (|47|) iff 
u = Q + v + Q~w and v, w are solutions of (fl5|) and (fl6|) respectively. Now 
using 1)44(1 we finish the proof. ■ 

Thus solution of Q with a of the form ([43(1 reduces to solution of the 
Schrodinger equations with quaternionic potentials f(|45jl and (|46|) which can 
be rewritten in a more explicit form as follows 

(-A- (a 2 -iDa x )I)v = (48) 

and 

(-A - ("a 2 + iDai)I)w = 0. (49) 
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Remark 13 In the case when a% does not depend on x\: a\ = a\{x2,x^) it 
is easy to see that equations \4&\ ) an d \4^i) are n °t independent. We have that 
if v is a solution of \4$ then w = ie\v is a solution of \4^i) an d vice versa. 
Using this fact we obtain a one-to-one correspondence between solutions of 
\4'l\j and solutions of i4&\ ). Such a correspondence is given by the operator 

n = i(/ + ie 1 -C + ie 1 C). 

It can be verified immediately that II 2 = / and that u is a solution of \47\) 
if and only if v = TLu is a solution of \4$ an d vice versa. That is when 
"a = ai(x2, xs)ei + 02^2 + 03^3 where 02, 03 are constants, solution of £7Jj 
reduces to solution of one Schrddinger equation with quaternionic potential 

Remark 14 Some classes of solutions of 14 £\ ) and \49[) can be obtained 
reducing the equations to scalar Schrddinger 's equations using the follow- 
ing idea (proposed in [15] in another setting). Consider equation \4$[) . If 
(iDoti — ~a 2 ) £ & or Dai £ © then we can look for solutions of \4<$ in 
the form v = (iDa\ + ~a 2 )f or v = (Dcti)f respectively where f is an un- 
known function. In the first case the equation reduces to the Laplace equation 
Av = and in the second to the Schrddinger equation (A + 7x 2 )v = 0. 

Suppose that neither of these is the case. Then denote (3 = iDoi\ and 
introduce j3$ as a scalar square root of [3 2 . The complex quaternions (3 = 
(3q + (3 and (3 = (3$ — (3 are conjugate zero divisors. Equation \4$\ ) then can 
be rewritten as follows 

(-A-((3 + ^ 2 )+f3)v = 0. 

Looking for its solutions of the form v = (3f we reduce it to the scalar 
Schrddinger equation (A + {(3q + He 2 ))v = 0. 

5 Conclusions 

We have shown that an ample class of physically meaningful partial differ- 
ential systems of first order are equivalent to a single quaternionic equation 
which in its turn reduces in general to a Schrddinger equation with quater- 
nionic potential, and in some situations considered in Subsection 4.1 this 
last can be diagonalized. The rich variety of methods developed for differ- 
ent problems corresponding to the Schrddinger equation can be applied to 
the systems considered in the present work. 
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